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Abstract 

We derive the optimal investment decision in a project where both demand and investment 
costs are stochastic processes, eventually subject to shocks. We extend the approach used in 
[T], chapter 6.5, to deal with two sources of uncertainty, but assuming that the underlying 
processes are no longer geometric Brownian diffusions but rather jump diffusion processes. 
For the class of isoelastic functions that we address in this paper, it is still possible to derive 
a closed expression for the value of the firm. We prove formally that the result we get is 
indeed the solution of the optimization problem. Finally, we derived comparative statistics 
for the investment threshold with respect to the relevant parameters. 


1 Introduction 

A real option is the right, but not the obligation, to undertake certain business initiatives, such 
as deferring, abandoning, expanding, staging or contracting a capital investment project. Real 
options have three characteristics: the decision may be postponed; there is uncertainty about 
future rewards; the investment cost is irreversible. There are several types of real options and 
different kinds of questions, but in all of them there is the issue regarding the optimal timing to 
undertake the decision (for instance, investment, abandon or switching). Consult [1] for a survey 
on real options analysis. 

In this paper we propose a two-fold contribution to the state of the art of real options analysis. 
On one hand we assume that the uncertainty involving the decision depends on two stochastic 
processes: the demand and the investment cost. On the other hand we allow the dynamics of these 
processes to be a combination of a diffusion part (with continuous sample path) with a jump part 
(introducing discontinuities in the sample paths). So, the model presented here expands the usual 
setting of options and adds new features, including the ability to model jumps in the demand 
process, and including stochastic investment costs. 
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Figure 1: Global Apple iPhone sales from 


quarter 2007 to 4*^ quarter 2014 (in million units) 


In the early works in hnance, regarding both hnance options and investment decisions, the 
underlying stochastic process modeling the uncertainty was frequently assumed to be a continuous 
process, as a geometric Brownian motion (GBM, for short). There are many reasons for this, but 
this is largely due to the fact that the normal distribution, as well as the continuous-time process 
it generates have nice analytic properties [2]. 

Observations from real situations show that this assumption does not hold. The bar chart in 
Figure shows the global sales numbers of iPhones from the 3'^'^ quarter of 2007 to the end of 
2014, with data aggregated by quarters (source: Statista 2014). It is clear that in this case upward 
jumps occur with a certain frequency, that may be related with the launching of new versions of 
the iPhone. 

Nowadays, with the world global market, there are exogenous events that may lead to a sudden 
increase or decrease in the demand for a certain product. Glearly, these shocks will considerably 
affect the decisions of the investors. Also the investment costs depend highly on other factors, 
such as the oil, raw materials or labour prices. Therefore, when we assume the investment to be 
constant and known, we are simplifying the model. This problem is even more important when 
one considers large investments that take several years to be completed (e.g. construction of dams, 
nuclear power plants, or high speed railways). 

The example regarding the number of sales of iPhones shows the existence of jumps in the 
demand process. In fact, empirical analysis shows there is a need to contemplate models that 
can take into account discontinuous behavior for a quite wide range of situations. For example, 
regarding investment decisions, in particular in the case of high-technology companies - where 
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the fast advances in technology continne to challenge managerial decision-makers in terms of the 
timing of investments [7] - there is a strong need to consider models that take into acconnt sndden 
decreases or increases in the demand level |3]. 

Morever, examples regarding jnmps in the investment are also easy to hnd. The jnmps may 
represent nncertainties abont the arrival and impact (on the nnderlying investment) of new in¬ 
formation concerning technological innovation, competition, political risk, regnlatory effects and 
other sonrces 13 E]. It is particularly important to consider jumps process in investment problems 
related with R&D investments, as |1] suggest. Furthermore, when one is deciding about invest¬ 
ments in energy, which are long-term investments, one should accommodate the impact of climate 
change policies effect on carbon prices, as [D] points out in their report. For instance, the Horizon 
2020 Energy Challenge is designed to support the transition to a secure, clean and efficient energy 
system for Europe. So, if a company decides to undertake any kind of such investment during 
this time period, probably will apply to special funds and thus the investment cost will likely be 
smaller. After 2020, there are not yet rules about the European funding of such projects. There¬ 
fore, not only the investment cost is random, but also will suffer a jump due to this new regulatory 
regime (see http://ec.europa.eu/easme/en/energy for more details). 

Following this reasoning, in this paper we assume that the decision regarding the investment 
depends on two stochastic factors: one related with the demand process, that we will denote 
by X = {Xt,t > 0}, and another one related with the investment cost, I = {It,t > 0}. We 
allow that both processes may have discontinuous paths, due to the possibility of jumps driven by 
external Poisson processes. Later we will introduce in more detail the mathematical model and 
assumptions. 

We also claim that we extend the approach provided in [Tj, in Chapter 6.5. In this book, the au¬ 
thors consider that there are two sources of uncertainty (the revenue flow and the investment cost), 
both following a GBM, with (possible) correlated Brownian motions. Then the authors propose a 
change of variable that will turn this problem with two variables into a problem with one variable, 
analytically solvable. Here we follow the same approach, but with two different features: hrst the 
dynamics of the processes involved are no longer a GBM and, secondly, we prove analytically that 
the solution we get is in fact exactly the solution of the original problem with two dimensions. In 
order to prove this result we derive explicitly the quasi-variational inequalities emerging from the 
principle of dynamic programming, obtaining the corresponding Hamilton-Jacobi-Bellman (HJB) 
equation. Moreover, using a verihcation theorem, we prove that the proposed solution is indeed 
the solution of the original problem [1316]. 

The remainder of the paper is organized as follows. In Section 2 we present the framework 
rationale and the valuation model, and in Section 3 we formalize the problem as an optimal 
stopping one, we present our approach and provide the optimal solution. In Section 4 we derive 
results regarding the behavior of the investment threshold as a function of the most relevant 
parameters and, hnally. Section 5 concludes the work, presenting also some recommendations for 
future extensions. 
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2 The model 


In this section we introduce the mathematical model along with the assumptions that we use 
in order to derive the investment policy. For that, we start by presenting the dynamics of the 
stochastic processes involved, namely, the demand and the investment expenditures processes. 
We also introduce the valuation model we use to derive the optimal investment problem that we 
address in this paper. 


2.1 The processes dynamics 


We assume that there are two sources of uncertainty, that we will denote by X and /, representing 
the demand and the investment, respectively. Moreover, we consider that both have continuous 
sample paths almost everywhere but occasionally jumps may occur, causing a discontinuity in the 
path. 

Generally speaking, let Y = {Yt : t > 0} denote a jump-diffusion process, which has two types 
of changes: the ’’normal” vibration, which is represented by a GBM, with continuous sample 
paths, and the ’’abnormal”, modeled by a jump process, which introduce discontinuities. 

The jumps occur in random times and are driven by an homogeneous Poisson process, hereby 
denoted by {N^, t > 0}, with rate A. Moreover, the percentages of the jumps are themselves random 

variables, that we denote by although they may be of different magnitude, 

they all obey the same probability law and are independent from each others. We let Ui > — 1, 
allowing for positive and for negative jumps. 


Thus we may express Yt as follows: 


Nt 

Ilfi + t'i). (1) 

where yo is the initial value of the process, i.e., Yq = yo, /i and a > 0 represent, respectively, 
the drift and the volatility of the continuous part of the process and {Wt,t > 0} is a Brownian 
motion process. Furthermore, the processes {Wt,t > 0} and {Nt,t > 0} are independent and are 
also both independent of the sequence of random variables We use the convention that if 

Nt = 0 for some t > 0, then the product in equation ([^ is equal to 1. Note that considering the 
parameters associated with the jumps equal to zero, we obtain the usual GBM. 


Yt = yo exp 


a 


—-— Am ) t + aWt 


As both the Brownian motion and the Poisson processes are Markovian, it follows that Y is 
also Markovian. It is also stationary and its order moments are provided in the next lemma 
(see Appendix (6.1) for the proof). 


Lemma 1. For Y = {Yt : t > 0}, with Yt defined as in ([^, f, s > 0 and k eN, 


IE [I'i.l Y,] = YT exp 


a 


y,+ {k — 1)— — Am ) fc -|- A (IE (1 -|- U) 


- 1 


where IE [.| Tj] denotes the conditional expectation given the value of the process Y at time t. 
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Returning to the investment problem: in the sequel we assume that both processes X and I 
follow jump-diffusion processes, such that 


and 


Xt = xo exp 


It = it) exp 


a 


Hx -^ — ^xmx ) t + (JxWt 


X 


n ( 1 +un 


2=1 


(Tt 


- t + (^I^t 


TV/ 


n ( 1 + ^ 1 ) 


2 = 1 


We use an index X or / for each parameter to distinguish between the two processes. Finally, we 
assume that the both processes are independent from each other. 

Next we introduce the valuation model, that dehnes the production function for the problem 
that we propose to address. 

2.2 The valuation model 


In this paper we propose to solve the question of when to invest in a new project, when both the 
demand and the investment are stochastic and subject to shocks. Furthermore, we assume that 
the investment may take some time to really starts to pay back, and therefore there may be a gap 
of n time units between the decision and the effectiveness of the investment. 

We also assume that the proht of the hrm depends on the stochastic process of the demand, 
and we use the following notation: Vq (W) is the profit of the hrm before (after) the investment. 
We consider an isoelastic function: 

for j G {0,1}, with ki > kq > 0 and 6 > 0. So the investment in the new project only changes 
the model coefficients, kq and ki, and not the elasticity parameter, 6. Furthermore, 

D{Xt) = {k)-ko)X^ (2) 

is the proht’s diherence before and after the investment. 

Denoting the risk-neutral discount rate by p, then the value of the hrm if the investment is 
decided at time r is given by: 

pT-\-n r-\-oo 

/ VoiXt)e-'^^dt+ / V){Xt)e-^Ut - Ire-P\ 

Jo J r+n 

as long as r < -fcx). 

The goal of the hrm is to maximize its expected value, with respect to the time to decide to 
invest, that is, the hrm seeks to maximize the following expected value 

pT-\-n ( /‘H-oo 


IE 

L^o 

_ -[pi(xo,io) 


Vo{Xt)e-PUt + 


V){Xt)e-P^dt-Ire-^^\x{r<+oo} 


' r+n 


(Xo — Xo, /o — *o) 


T*T+n 


Vo{Xt)e-PUt + 


L^o 


'*+00 


V){Xt)e-PUt - Ire-P^ X{ r<+oo} 


' r+n 
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where in the rest of the paper denotes the expectation conditional on (Xq, Jo) = {xo,io)] 

moreover X{A} represents the indicator function of the set A and r is a stopping time for the 
hltration generated by the bidimensional process (X, I). 

Using (|^ and simple manipulations, we can re-write the previous expectation as 


]g(2:0,*o) 


b+oo 


Vo{Xt)e-P^dt 


Uo 


_|_ jg(a:o,io) 


'■+00 


D{X,)e-<^^dt - Ire-P^ X{ r<+oo} 


• r+n 


As the hrst part of the expression does not depend on the time r to invest, then the goal is to 
maximize the following functional: 




^+00 
' T+n 


D{Xt)e-P^dt-Ire-^^\x{ r<+oo} 


which represents the value of the option when the hrm exercises it at time r, with r < -t- 00 , given 
that the initial state is (xo,io). Following [^, we call J'^ the performance criterion. 


Note that in order to have r < -fcxD, one needs to assume that Vo{Xt)e ^^dt 

-l-cxD, which happens if and only if 

h > 0, 

where 


< 


a 


h — p — [ px + (d — 1)-^- ^x'^x ] 0 — Xx 


IE 


{1 + UxY 


- 1 


(3) 


since 


]gUo,*o) 


b+oo 


Vo{Xt)e-P^dt 


Uo 


^+00 


= KqXq 


—ht 


dt. 


From now on, we will be assuming such restriction in the parameters, to avoid trivialities. Next 
we state and prove a result related with the strong Markov property of the involved processes, 
that we will use subsequently in the optimization problem. 

Proposition 1. The performance criterion can he re-written as follows: 

r{x,z) = [e-^^g{Xr,Ir)X{r<+o.}] , 

with 


g{x, i) = {ki — kq) Ax^ — i 


and 


^—hn 


A = 


h 


( 4 ) 


Proof. Applying a change of variable in the performance criterion, we obtain 

/*+O0 


X{x,i) = 


:,-P^ 


-Ir>X{ T<-|-00} 
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Then, using conditional expectation and the strong Markov property, we get 






IE 


^ + 00 




{XrJr) 


r X{r<+oo} 


By Fubini’s theorem and the independence between X and /, it follows that: 


IE 


"+00 




L-'o 


{X^Jr) 


<•+00 


M[D{X^+t+n)\Xr\ 


From Lemma and from simple calculations, we have: 

IE [D{Xr+t^n)\Xr\ = («1 - Ko) Xl 
Plugging ([^ on (|^ results in: 


IE 


b+oo 


D{Xr+t+n)e-P^^^^^dt 


IJo 


{XrJr) 


— (Ki — Kq) X^ 


-,—hn 


h 


Therefore, in view of the dehnition of J'^ and A, we get the pretended result. 


(5) 

( 6 ) 


□ 


3 The optimal stopping time problem 


The goal of the hrm is to hnd the optimal time to invest in this new project, i.e., for every a: > 0 
and i > 0, we want to hnd V{x, i) and t* & T such that: 

V{x, i) = sup J^{x, i) = (x, i), {x, i) G M"*" x M"*" 

r&T 

where T is the set of all stopping times adapted to the hltration generated by the bidimensional 
stochastic process (X, I). The function V is called the value function, which represents the optimal 
value of the hrm, and the stopping time r* is called an optimal stopping time. 

Using standard calculations from optimal stopping theory (see m), we derive the following 
HJB equation for this problem: 

min{pV{x,i) — CV{x,i),V{x,i) — g{x,i)} = 0, V (x, f) G M’*'x M’*', (7) 


where C denotes the inhnitesimal generator of the process (X, I), which in view of [6], it is equal 
to 


CV{x,i) = 


d‘^V{x,i) 


O'j .0 


5W( 


x,l] 


+ (hx - X 


dV{x,i) 


+ {pi - Ximj) i 


.dV{x,i) 


dx'^ 2 di"^ ' ^ • \r~i 11 ) 

+Xx [lEt;x (U(a;(l + U^),i)) - V{x, ^)] + A, [lE^. {V{x, ^(1 + U^))) - V{x, 0] 


Note that in this equation we used the notation lE^x and lEjy/ to emphasize the meaning of such 
expected values; for instance Wijji means that we are computing the expected value with respect 
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to the random variable . Whenever the meaning is clear from the context, we will simplify the 
notation and will use just IE instead. 

The HJB equation represents the two possible decisions that the firm has. The hrst term of the 
equation corresponds to the decision to continue with the actual project and therefore postpone 
the investment decision. The second term corresponds to the decision to invest and hence the end 
of the decision problem. Thus 

C = {(x,i) e M’*' X : V{x,i) > g{x,i)] 

is the set of values of demand and investment cost such it is optimal to continue with the actual 
project and, for that reason, is usually known as continuation region, whereas in the stopping 
region, S, the value of the firm is exactly the gain if one decides to exercise the option. The sets 
C and S are complementary on x M’*'. We define 

r* = inf {f > 0 : (W, h) i C} , 

which is the hrst time the state variables corresponding to the demand and the investment cost 
are outside the continuation region. Note that Q means that 

pV(x, i) — CV{x,i) >0 A V{x, i) > g{x, i), V {x, i) G R'^ x M+. 


Moreover, 


V{x, i) > g{x, i) A pV(x, i) — CV[x, i) = 0, V [x, i) G C 


whereas 


V[x, i) = g{x, i) A pg{x, i) — Cg{x, f) > 0, V {x, i) G S. 

The solution of the HJB equation, V, must satisfy the following initial condition. 


i/(0+,i)=0, ViGM+, (8) 

which rehects the fact that the value of the hrm will be zero if the demand is also zero. Furthermore, 
the following value-matching and smooth-pasting conditions should also hold 


V{x,i) = g{x,i) and W{x,i) = Vg{x,i) 


(9) 


for {x, i) G dC, where dC denotes the boundary of C, which we call critical boundary, and V is the 
gradient function. In this particular case, with two sources of uncertainty, the critical boundary 
is a threshold curve separating the two regions. 


So in order to solve the investment problem, we need to hnd the solution to the following 
partial differential equation: 


pV (x, i) — CV (x, i) = 0 


( 10 ) 


in the continuation region and, simultaneously, to identify the continuation C and stopping S 
regions. 


Both questions are challenging: on one hand the solution to (10) cannot be found explicitly, 
and on the other hand we need to guess the form of the continuation region and prove that the 
solution proposed for the HJB equation satishes the verihcation theorem. Hereafter, we propose 
an alternative way to solve this problem, that circumvent the difficulties posed by the fact that 
we have two sources of uncertainty. 



3.1 Splitting the state space 

In this section we start by guessing the shape of the continuation region and latter we will prove 
indeed that our guess is the correct one, using the verification theorem. 

Our guess about the continuation region comes from intuitive arguments: high demand and 
low investment cost should encourage the hrm to invest, whereas in case of low demand and high 
investment cost it should be optimal to continue, i.e., postpone the investment decision. See Figure 
for the general plot of both regions. 


i 


C 


S 


-► 

X 

Figure 2: Space of process (X,/), C (continuation region) and S (stopping region). 


But we need to define precisely the boundary of C and for that we use the conditions derived from 
the HJB equation. We consider the set 


f/ = {(a:, i) G M’*' x : pg{x^ i) — Cg{x, i) < O} 


which, taking into account the expression of C (the inhnitesimal generator) and the g function, 
can be equivalently expressed (after some simple calculations) as 


U 


|(a;, i) e x 


x^ p — Pi 

i (ki — Ko) Ah 


By Propositions 3.3 and 3.4 of [B], we know that U O C. Further, if 17 = 0 then V{x,i) = 
g{x,i), V {x,i) G x M"*" and r* = 0. This case would mean that it is optimal to invest right 
away. So we need to check under which conditions this trivial situation does not hold. We have 
(K,^-KQ)Ah ^ 0 if only if p — p/ < 0, as the denominator is positive; as ^ > 0, this would imply 
that 17 = 0 and r* = 0. Therefore, if p/ > p the firm should invest right away, which is coherent 
with a hnancial interpretation: if the drift of the investment cost is higher than the discount rate, 
the rational decision is to invest as soon as possible. To avoid such case, we assume that this does 
not hold, i.e., 

p > pi. (11) 
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Next we guess that the boundaries of U and C have the same shape, and therefore C may be 
written as follows: 


C 


|(a;,i) e M+ X M+ 


: 0 < 



5 


where g* is a trigger value that satishes 


q* > 


P- 

(ki — Ko) Ah 


( 12 ) 


In Figure the thick line represents the boundary of the set U and the other dashed lines are 
some possible boundaries of the continuation region, C (as U which would be the case). 



Figure 3: Set U (bold line) and some possible boundaries of the continuation region, C. 


3.2 Reducing to one source of uncertainty 


It follows from the last section that the decision between continuing and stopping depends on the 
demand level and investment cost only through a function of the two of them: 


X 


e 



Then the g function present in the performance criteria can also be written in terms of this new 
variable, as g{x,i) = with 

liq) = {ki - kq) Aq - 1 (13) 


and, in view of (13), we propose that 


V{x,i) =if {q), 


where / is an appropriate function, that we need to derive still. 

In chapter 6.5 of [1] is discussed a problem with two sources of uncertainty following a GBM. 
Using economical arguments, they propose a transformation of the two state variables into a 
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one state variable, and thus reducing the problem to one dimension, that can be solved using a 
standard approach. 

In the problem that we address we don’t use any kind of economical arguments but we rely on 
the definition of the set U. Furthermore, the state variables that we consider in this paper follow 
a jump-diffusion process, and thus the dynamics is more general than the one proposed by [T]. 
Therefore, we may see our result as a generalization of the results provided by pQ. 

Following the approach proposed in the previous section, we write the original HJB equation 
now for the case of just one variable: 


min {pf{q) - CQf{q), f{q) - l{q)} = 0, V g G 


(14) 


where the inhnitesimal generator of the unidimensional process Q = 


Q = {Qt 


X 


= \ is 


+ "{(!) + 


a 


Px + (0 — 1 )-^- i^xnix ] 9 — {pi — Ximi) 


+ [ii^i ~ Ximj) — {Xx + A/)] /(g) 


-|-AxlE 


f(q{l + U^y 


A/IE 




{i + u‘)fU{i + u‘) 


‘if'ii) 

(15) 


where this equation comes from C definition, the relationship between V and /, and simple but 
tedious computations. 

The corresponding continuation and stopping regions (hereby denoted by Cq and Sq, respec¬ 
tively) can now be written as depending only on g: 

^<9 = ^ : 0 < ? < ?*} and = {g G : g > g*} , 

and thus g* is the boundary between the continuation and the stopping regions, as it is usually the 
case in a problem with one state variable. We note that now we are in a standard case problem 
of investment with just one state variable (see, for instance, chapter 5 of P) but with a different 
dynamics than the usual one (the process Q is not a GBM). 

Moreover, the following conditions should hold 


/(0+)=0, /(g*) = /(g*) and f{q) = l\q)\q=g*- (16) 

Next, using derivations familiar to the standard case, we are able to derive the analytical solution 
to equation (14). 


Proposition 2. The solution of the HJB Equation (14) verifying the conditions (16), hereby 
denoted by f, is given by 


fiq) = ^ 


ro 


0 < g < g* 
{ki - Ko)Aq-l q>q* 


where 


q* = 


ro 


{ki - kq) a (ro - 1 ) 


(17) 


(18) 
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and tq is the positive root of the function j : 


J[r) = 


ax „2 
2 2 


^2 _|_ 


Hx -^ — ^xmx j ^ ~ ^ ~ Ximi 


+ ifij - Xjmj) - (Ax + Xj)-p + AxE [(1 + U^Y^] + A,E [(1 + • (19) 

Proof. For the case q > q* (corresponding to the stopping region), the value function is equal to I, 


the investment cost, dehned on (13), and thus the result is trivially proved. We then must prove 


the result for the continuation region, for which q < q*. 

In this case, we propose that the solution of the HJB in the continuation region is of the form 


C(g) = bq'^°, where b and ro are to be derived. As C(0) = 0 (from the condition (16)), it follows 
that ro needs to be positive. 


Furthermore, in view of the dehnition of Cq (15), if we apply it to the proposed function it 
follows that: 

pC{q) - = -bq''°j{ro), 


where j is dehned on (19). Therefore, in the continuation region 

pCiq) - ^gCiq) = 0 j(ro) = 0 , 

meaning that tq is a positive root of j. Next we show that j has one and only one positive root 
and thus tq is unique. The second order derivative of j is equal to 


j"(r) = a^e^ + a] + XxO^lE [ln(l + U^)] ^ (1 + U^Y^ 


+ A/E 


[\n{l + UYY {1 + UY^- 


Given that the volatilities and the elasticity are positive, > —1 a.s. and > —1 a.s., it follows 
that j"{r) > 0, Vr G M. Thus, j is a strictly convex function. Further, limr^+oojY) — 


j(0) = Pi — p < 0 (by assumption (11)), which means that j(0) < 0. Then, since j is continuous, 
it has an unique positive root. 


It remains to derive b. For that we use conditions presented in (16), and straightforward 
calculations lead to: 

1 / 1 ^ 
b = - 


ro 




where q* must be given by (18). In order to ensure that this leads to an admissible solution for q* 


(as q* must be positive), one must have ro > 1. 

The remainder of the proof is to check that 0 is the solution of the HJB Equation ( [l^ . For 
that we need 

i) To prove that in the continuation region f{q) > l{q). In the continuation region f{q) = 

^ (-X ) • So, we dehne the function 

0-1 \q* J ’ 


ro- 


p{q) = 


q"" - {ki - Ko) Aq + I 


ro - 1 \q 
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with derivatives 


^'{q) = ^ - («i - Ko) A and (p"(q) = Tq ^ 

Since (fi"(g) > 0, is a strictly convex function. As (fi'(g) = 0 g = g*, then g* is the 
unique minimum of the function ip with p (g*) = 0. 

Summarizing, the function p is continuous and strictly convex, has an unique minimum at 
g* and p (g*) = 0, then p{g) > 0, Vg G M’*'. Therefore, /(g) > /(g), for Vg G M. 

ii) To prove that in the stopping region p/(g) — Cqf^g) > 0. 

In the stopping region /(g) = /(g), then 


pf{q) - ^Qfiq) = - «o) Ahg + {pi - p). 


Given that g > g* and using the inequality (12), we conclude that 

P- Pi 


pC{g) - ^gCiq) > («i - «o) Ahg* + {pj - p) > (ki - Kq) Ah 


{ki — Kq) Ah 


+ Pj - p = 0. 


Therefore, we conclude that the function / is indeed the solution of the HJB Equation (14). □ 


We have already seen that in order to avoid the trivial case of investment right away, one must 
have p > Pi. Moreover, from the previous proof of the proposition, one must also assume that 
vq > 1, which implies that j(l) < 0, condition equivalent to h > aj, where h is dehned in (|^. 

Assumption 3.1. ITe assume that the following conditions on the parameters hold: 

p > Pi and h > aj. 


3.3 Original problem 


In this section we prove that the solution of the original problem can be obtained from the solution 
of the modihed one. 


Corollary 1. The solution of the HJB Equation ([^ verifying the initial condition ([^ and the 
boundary conditions (|^, hereby denoted by V, is given by 


V{x,i) = 


(ki-ko)Ax'^ Y° ^ rp-l 


ro 


(ki — Ko) Ax^ — i 


o<Y <q* 

^>q* 


( 20 ) 


for (x,i) G 1^^ X IR''", where q* is defined on (18) and rg is the positive root of the function j, 
defined on (19). 
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Proof. Considering the relations of the functions g and / with the function V, it follows that the 
function V proposed in Equation (20) is the obvious candidate to be the solution of the HJB 
Equation ([^. So next we prove that indeed this is the case. In the following we use some results 
already presented in proof of Proposition 


Let {x,i) G x M+ and q = ^. Using some basic calculations one proves that CQV{x,i) = 
iCgf^q) and pV{x,i) — CV{x,i) = i [pf{q) — Cgfi^q)]. Since i is positive, the signal of pV{x,i) — 
CV{x,i) is exactly the signal of pf{q) — CQf{q). By construction, pf{q) — Cqfi^q) = 0 in the 
continuation region and pf{q) — C,f{q) > 0 in the stopping region. Therefore, we just need to 
prove that in the continuation region V{x,i) — g{x,i) > 0. This follows immediately because 
V{x,i) — g{x,i) = i [f{q) — i{q)], i is positive and we previously proved that f{q) > l{q). 


Finally, from the initial and boundary conditions of the transformed problem (16), U(0’'',i) = 
= 0 and for {x*,i*) G dC, i.e., ^ = q*, we obtain the boundary conditions of the original 
problem (|^. □ 


4 Comparative statics 


In this section we study the behavior of the investment threshold as a function of the most relevant 
parameters, in particular the volatilities ax and ai, the intensity rates of the jumps. Ax and A/, 
and the magnitude of the jumps, mx and m/. In order to have analytical results, we assume that 
the magnitude of the jumps (both in the demand as in the investment) are deterministic and equal 
to mx and m/, respectiveljj^ 

Before we proceed with the mathematical derivations, we comment on the meaning of having 
g* monotonic with respect to a specihc parameter, as we are dealing with a boundary in a two 
dimensions space. For example, assuming that g* increases with some parameter y, this would 
mean that g*(2/i) < g*( 2 / 2 ), for yi < y 2 (leaving out the other parameters). Consider then the 
corresponding continuation regions: 


Cl = {x,i) G 


X 


X 


: 0 < — < g* 


and 


Ca = S {x,i) G 


X 


X 


: 0 < — < ga ^ 


with g^ = q*{yi) and g^ = g*(|/ 2 )- In Figure]^ we can observe the graphic representation of the 
boundaries which divide the stopping and the continuation regions, for each case, and, as expected, 
we would have Ci C Ca. It follows from the splitting of the space into continuation and stopping, 
that you in this illustration, we would stop earlier for smaller values of the parameter. 


^Numerical results suggest that the distribution of the jumps is not very relevant in terms of its effect on the 
investment thresholds. 
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Figure 4: The full (respectively, dashed) line is the boundary of the Ci (respectively, C 2 ) set. 

In the previous sections we consider j only as a function of r. As in this section we study 
the behavior of j with respect to the parameters of both processes involved, we change slightly 
the notation, in order to emphasize the dependency of j in such parameters. Therefore, we 
define the vectors (px = 4>i = (h/) ct/. A/, m/) , = {9,p,<px) and finally 

Q = {0,p, <px, <Pi), that we use in the following dehnition: 


jir,Q) = 


a 




CX7 


+ 


9'X 


— Xx'^x \ P'l — w ~ Ximj 




+ {pi — Ximi) — (Ax + Xi) — p + Ax(l + rnx) ^ + A/(l + mj) 


\ 1 — r 


Furthermore, we let ro(^) denote the positive root of j and we change accordingly the dehnition 
of h as follows: 


his) = p - i^px + {0 -l)^ - Axmx^ 6^ - Ax (^(1 + mxf - l) 


and hnally 


q\Q) = 


ro{Q) 


{ki - Ko) A(^) ro{Q) - 1’ 

with A(<j) dehned on (|^. 

We start by proving the following lemma: 

Lemma 4.1. The investment threshold changes with the parameters according to the following 
relation: 


dq*{Q) 

dy 


' Aro(g)(ro(e)-l)(nfe(^)+l) dhp) , dj{r,Q) 
h{<;) dy dy 


r=ro 


A{ki - Ko)A{^){ro{Q) - 1)2 


X 


dj{r,Q) 

dy 


ify e (px 


ify e cpi 


r=ro 


with A = 

dr 


Furthermore, as A > 0, then the sign of the derivative of q* depends only 


r='ro 


on the derivatives of h and j with respect to the study parameter. 


15 





















Proof. As 


dq*{Q) 

dy 


1 

(ki - Ko)A(<j)(ro(£>) - 1)2 


?'o(g)(''o(g)-l)(’^fe(‘;)+l) dh(i;) 

h{<;) dy 

dro{Q) 

dy 


dro(Q) 

dy 


iiy ^ (f>x 
if 1/ e 07 


the result follows from the implicit derivative relation 


dro(Q) 

dy 


dj(r,g) 

dy 

9j{r,Q) 

dr 


r=ro 


It remains to prove that A > 0. As a function of r only, j is continuous and strictly convex, 
with limr^^ooj(r) = + cxd , j(0) < 0 and j(l) < 0, as we have already stated. Since tq is its positive 
root, in view of the properties of j, it follows that it must be increasing in a neighborhood of tq, 
and thus A > 0. □ 


In the following sections, we study the sign of each derivative, whenever possible. Our goal is 
to inspect if the trigger value, q*, increases or decreases with each parameter in (px and 0/. 


4.1 Comparative statics for the parameters regarding the demand 


In this section we state and prove the results concerning the behavior of q* as a function of the 
demand parameters. Before the main proposition of this subsection, we present and prove a 
Lemma that provide us the admissible values for the parameters. 


Lemma 4.2. In view of the Assumption 5A 
hold. The drift is upper bounded, i.e., 


the following restrictions on the demand parameters 




a 


P- ( {0 -1)— -\xmx 


^ 6^ - Ax ((1 + mx)^ - a] 


The domain of the volatility parameter depends on 6. If 0 < 6 < 1, ax is lower hounded, ax > B, 
otherwise is upper hounded, < ax < B, where 


B = \ max < 0, 


6(6-1) 


P - (px - Ax^x) 6 - Xx [{1 + mx) -l) - a] 


Relatively to the rate of the jumps, the domain also depends on 6. If D < 6 < 1, Xx is lower 
hounded, Xx > C, otherwise is upper bounded, D < Xx < C, where 


P- {px +(0 -1)^^6 - a]^ 

(l + mxY - (6mx+ 1) j 


Regarding the jump sizes, we have the following: let and ^2 be the negative and positive solutions 
of the equation 


(1 + xY — (6x + 1) 


p - {px + (6^ - 1) ^ j 6 -aj 
Xx 


( 21 ) 
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respectively, when there are solutions. Moreover, let 


For 0 < 9 < 1, 


D 


p- (^Fx + {9- 1)^) 9-aj 
Ax 


• if D < 6 — 1 then mx > ^ 2 / 

• ifO — l<D<t) then mx G (—1, .^i) U (.^ 2 , + 00 ); 

• if D > 0 then mx > — 1 . 

For 6 > 1, 

• ifO<D<9 — 1 then mx G (■Ci; 'C 2 ); 

• if D > 9 — 1 then mx G (—1, ^ 2 ) [5 


and if 9 = 1 and D > 0, then mx > —1- 


( 22 ) 


Proof. The sets of admissible values for the parameters px and ax come from straightforward 
manipulation of the second restriction of Assumption 3.1[ 

The results regarding the jumps parameters (Ax and mx) rely on the properties of the function 


T„(a;) = T„,2(a;) - To,i(a;), (23) 

with To i(a;) = ax + 1 and Tap^x) = (1 + x)", with x > — 1 and a > 0. The graph of the 
function is a tangent line of the graph of the function To^2 on the point (0,1)- Also, as 
^a2(^) = ~ 1)(1 + the function Tq,_ 2 is convex if a > 1 and it is concave otherwise. 

Then, for x > — 1, Tapix) > To,i(x) if a > 1 and To_i(x) > Ta,2{x) if 0 < a < 1. Then it follows 
that Tq,(x) > 0 if a > 1 and To(x) < 0 if 0 < a < 1. 

With respect to Ax, in view of the restriction h > aj, 

AxTe(mx) < P — (^Px + (9 — 1)“^^ 9 — aj 

and therefore the restriction on this parameter follows. 

Proceeding with the analysis with respect to the jump size, mx, we start by noting that 
T.(-l) = a — 1, Tq( 0) = 0 and lima;^oo Tq,(x) = +C)oy;{Q,>i} — ooX{a<i}. From these properties 
we know the shape of the function Tq, for each a. The restriction h > aj is equivalent to have 
Tg{mx) < F) and therefore standard study of functions lead us to the results presented in the 
Lemma. 


^The case D < 0 is equivalent to h < aj, 


and therefore Assumption 


3.1 


does not hold. 


□ 
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Now we are in position to state the main properties of q* as a function of the demand param¬ 


eters. From Lemma 4.1, we can only get results about the behavior of q* as a function of the 


parameters if 


dhji;) 

dy 


and 


dj{r,Q) 


dy 


have the same sign. Therefore we are not able to present a 


r=rQ 


comprehensive study of q* for all values of 9. 

Proposition 3. For F 9 < 1, the investment threshold q*: 

• increases with the demand volatility, ax and jump intensity, \x; 

• is not monotonic with respect to the magnitude of the jumps in the demand, mx ■ In fact, 

~ if D < 9 — 1 , then q* increases for mx > ^ 2 ! 

— if 9 — 1 < D <0, then q* decreases for mx € (—l,^i) and increases for mx > ^ 2 ,' 

— if D > 0, then q* decreases for mx G (—1, 0) and increases for mx > 0, 


where D is defined on (22) and .^1 and ^,2 are,respectively, the negative and positive solutions 


of the Equation (21), when there are solutions. 


Proof. The proof for ax is straightforward knowing the derivatives 

= 6ro{Q){9ro{g) - l)ax. 


^ mr, e) 


da 


X 


da 


X 


r=ro 


Continuing with Ax, we calculate the derivatives 
dh{(; 


dXx 


dj{r, g) 


dX 


X 


= 9mx -Fl - (1 + mxY = -Tg{mx) 

= (1 + mxY''°^^^ - {9ro{g)mx + 1 ) = ^eroiQ){nrx), 


r=ro 


with Tq, given by (23). Taking into account the properties of such function, we conclude the result 
for Ax- 

Finally, for mx we also calculate the derivatives 
dh{(; 


dmx 


dj{r, g) 


dm 


X 


= 9Xx [1 — (1 -I- mx)^ = 9Xx^e{niix) 

= -dnQ{g)Xx [1 - (1 + = -9ro{g)Xx^eTo(Q){'dax) 


r=ro 


where 

$«(a;) = l-(l-Fxr, (24) 

with X > —1 and a > 0. As $^(2;) = — a(l -1- is decreasing if a > 0, constant and equal 

to zero if a = 0, and increasing if a < 0. Moreover, as $a(0) = 0, then: 
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• if a > 0, $Q, is positive when —1 < x < 0 and is negative for x > 0; 

• if a < 0, is negative when —1 < x < 0 and positive for x > 0. 

Given these results, the proof regarding mx is concluded. □ 


Remark 4.3. Numerical examples show us that the investment threshold q* is not monotonic with 
respect to the demand’s drift, qix- We cannot have analytical results for the behavior of q* with 
respect to jxx as the signs of the involved derivatives are opposite: 

0ro{Q). 


= —6 and 

OfiX 


dfix 


r=ro 


4.2 Comparative statics for the parameters regarding the investment 


Now we state and prove the results concerning the behavior of q* as a function of the investment 
parameters. In view of Lemma 4.1, we just need to assess the sign of the derivative of j with 


respect to each one of the parameters. 


Proposition 4. For values of investment parameters such that Assumption 3.1 holds, the invest¬ 
ment threshold q* 


• decreases with the investment drift pij; 

• increases with the investment volatility, aj and jumps intensity, Xj; 

• is not monotonic with respect to the magnitude of the jumps in the intensity; it decreases 
when mi e (—1,0) and increases when mi > 0. 


Proof. The proof for /x/ and cx/ is trivial. In order to pove the result for A/, we study the behavior 
of the function 0: 

Q(^) ^ 9j^, q) ^ _ X) + (1 + miY~'' - 1. 

uXj 

We compute the second derivative, 

0"(r) = [ln(l + miff (1 + > 0, 


whence 0 is a strictly convex function. Furthermore, 0(0) = 0(1) = 0, which lead us to conclude 
that 0(r) > 0, for r > 1. Therefore, = 0(?"o(^)) > 0, and the result follows. 


r=rQ 


Finally, in order to study the behavior of q* with respect to mi, we need to assess the sign of 
the following derivative: 

dj{r, q) 


dmj 


= A7(ro(p) - 1) [l - (1 + mi) . 

Taking into account the function defined on (24), we notice that 


r=ro 


dmj 


= A7(ro(p) 


r=rQ 


l)4>_ro(e)(^v)- As ro(p) > 1 and in view of the properties of function the result follows. □ 
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5 Conclusion 


In this paper we proposed a new way to derive analytically the solution to a decision problem 
with two uncertainties, both following jump-diffusion processes. This method relies on a change 
of variable that reduces the problem from a two dimensions problem to a one dimension, which 
we can solve analytically. 

Moreover, we also presented an extensive comparative statics. From such analysis, we have 
concluded that the behavior of the investment threshold with respect to demand and investment 
parameters is similar, i.e., the impact of changing one parameter in the demand process is qual¬ 
itatively the same as changing the same parameter in the investment process. This result holds 
for all except for the drift: the investment threshold is monotonic with respect to the investment 
drift but this does not hold for the demand drift. 

For some parameters, the result is as expected, and similar to the one derived in earlier works 
(like the behavior of the investment threshold with the volatility of both the demand and the 
investment), whereas in other cases the results are not so straightforward. For instance, the 
impact of the jump sizes of the demand depend on some analytical condition whose economical 
interpretation is far from being obvious. 

Our method may be extended for more general proht functions. In our research agenda we 
plan to study properties of the functions such that the approach proposed in this paper will still 
be valid and useful. 


6 Proofs 


6.1 Proof of Lemma [T] 

Proof. Applying expectation and re-writting Y^+s, we have 


IE[r4.|F,] =y,‘exp 


a 

Y 


] ks 

IE 

/ 



Nt 


t + S 


exp[ak{Wt+s-Wt)] n (1 + ^c 


i=Nt+l 




Using the fact that the Brownian motion and the Poisson process have independent and stationary 
increments, and that the processes are independent, we obtain 


IE [I'i.l y,] = n^exp 


a 


/i —-— Am ) ks 


IE [exp {akWs)] IE 


Ns 




2 = 1 


( 25 ) 


As U, then using the probability generating function of a Poisson distribution and 

the tower property of the conditional expectation, we conclude that 


IE 


Ns 


n (1+ Ut)‘ 


2=1 


= exp 


As IE 


(1 + f/)' 


- 1 


( 26 ) 
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The result follows using the moment generating function of the Normal distribution, and replacing 
® in d^. □ 
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